#ifndef MD_H_INCLUDED
#define MD_H_INCLUDED

double MD(particula *part, parametros param)
{
    int i, j;
    double vcm[3];
    double epoten;
    int cual=1;

    ///Calculamos las fuerzas
    epoten = interaccion(part, param);
    for(i=0;i<param.npart;i++)
    {
        ///Ecuaciones del movimiento
        for(j=0;j<3;j++)
        {
            part[i].vel[j]+=part[i].force[j]*param.Dt/part[i].mass;
            part[i].pos[j]+=part[i].vel[j]*param.Dt;
            boundary(&part[i],param);
        }
    }
    rescaleVel(part, param);
    return epoten;

}


#endif // MD_H_INCLUDED
